How big is the source that produces quark gluon plasma in heavy ion collisions? 
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PACS: 12.38.Mh, 25.75.-q,24.85.+p, 25.75.Gz 

We study, for the first time, the spatial extension of the 
"source" that produces quark gluon plasma (QGP) in ultra 
relativistic heavy ion collisions (URHIC). The longitudinal di- 
mension is studied as a function of time as the system evolves. 
The source size is found to exhibit a novel non-classical fea- 
ture. 



I. INTRODUCTION 

The ultimate aim of studying ultra relativistic heavy 
ion collisions (URHIC) is not merely to establish the pro- 
duction of Quark Gluon Plasma (QGP) - by looking at 
various signals - but to get a complete picture of the 
space-time evolution of QGP. In short, one hopes to de- 
scribe the evolution of this unique deconfined state of 
hadronic matter in all its stages - production, equilibra- 
tion and hadronisation. Ambitious that this might seem, 
it has indeed proved impossible to separate the study of 
signals such as J/^* suppression or strangeness enhance- 
ment virithout having an at least approximate idea of the 
production and evolution of QGP. 

We have recently proposed a mechanism [1] for produc- 
ing QGP in URHIC - the soft and the semisoft quarks 
and gluons that constitute the bulk of the plasma. The 
mechanism manifests as the source term in the transport 
equation for the distribution functions for the partons. 
This source term, which is the rate at which the partons 
are produced in the extended phase space, has the fol- 
lowing features: (i) it is defined in the two particle phase 
space, (ii) the production rate is non-Markovian in time 
[3], (iii) the vacuum has a dynamical role to play [2] - 
it acts both as a source and a sink, (iv) the partons are 
quasiparticles with a finite life time [4] , (v) the rate eval- 
uation takes into account the time scales involved in the 
production vis-a-vis the time intervals over which the rate 
is determined, (vi) the phase space dependence does not 
violate the uncertainty principle, and finally, (vii) the dy- 
namical nature of the colour charge is manifest. Given all 
these attributes, it is natural to ask whether the source 
term can throw light on the spatial dimensions of the 
QGP as it evolves in time. Although a complete answer 
to this question has to be kept pending until the trans- 
port equation is solved in a self consistent manner, the 
purport of this paper is to show that the source term by 
itself does have interesting information. We shall illus- 
trate this by looking at the longitudinal dimensions in a 
generic example. 



One promising observational tool for studying dynam- 
ically the spatial extension of the fireball is the measure- 
ment of HBT correlations [5-7] in two and three pions 
(and kaons) that are produced in URHIC. In particle 
physics, this idea was first proposed by Goldhaber et. al. 
[8] for pp collisions. In the context of QGP, a major im- 
petus for such a study is the hope that a knowledge of 
the system size would also settle whether the source pro- 
ducing the pions/kaons is nuclear, or its deconfined state. 
A useful tool that it is, measurements in HBT are still 
beset with many uncertainties - of extracting and inter- 
preting data. Much of it is due to the fact that the source 
here is not static; it evolves rapidly, over time scales of 
several fermis and over length scales of the same order. 
The reconstruction of the source history from HBT ob- 
servations is, therefore, not straight forward: it has been 
pointed out [9,10] that what gets determined is not the 
true " size" of the system, but only the so called regions of 
" space time homogeneity" . The other complications are 
the finite life time of the source, inhomogeneous tempera- 
ture profiles and a strong collective dynamical expansion 
[9]. At a more basic level, HBT studies have the best 
utility if the pions are produced from a chaotic source; 
there is no clinching justification (theoretically) to as- 
sume that in the case of URHIC; studies based on the 
Lund model [11,12] indicate that the sources do possess 
non-chaoticity. An empirical test to check the chaoticity 
is to measure the so called A parameter in two pion cor- 
relations. However, A itself depends on other paramters 
which are either incompletely or inaccurately known [13]. 

Apart from these caveats, the extraction of radii in- 
volves additional assumptions: the interpretation and 
parametrisation of HBT data assume a Gaussian pro- 
file for the source. Reasonable and natural that it seems, 
it still needs justification. Next, the (spurious) contri- 
butions from the Coulomb interaction have to be sub- 
tracted. More inportantly, one needs a criterion for 
"emission instant" and the "duration of emission of the 
particles" - both of which are ill defined in quantum 
mechanics [1,10]. The phase space coordinates should 
also be sufficiently coarsened so as not to violate uncer- 
tainty principle [1]. These problems are nontrivial to 
handle both in experiments and theoretical simulations 
[10]. An illustrative instance of the present situation is 
the so called RHIC puzzle [14] which is still not under- 
stood properly. 

In short, the utility of the HBT analyses gets enhanced 
if they can be supplemented with an independent theoret- 
ical investigation of the space time dynamics of QGP. As 
mentioned at the beginning, the transport studies pro- 
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vide the required frame work. We shall show that the 
the source term proposed in [1] does indeed give valuable 
information on the spatial extension of QGP. 

The next section introduces the source term derived 
in [1] and the subsequent section discusses the spatial 
extension of the source. For the sake of simplicity, we 
consider only the longitudinal dimension. 

II. THE SOURCE TERM 

We refer to [1] for the motivation and details of de- 
riving the source term. To summarise in a nut shell, 
the production of QGP takes place via the decay of a 
mean chromoelectric field (CEF) that is produced in be- 
tween the two nuclei soon after they have collided and 
start receding from each other. By energy momentum 
conservation, the CEF acquires a space - time depen- 
dence, thanks to which the instability of the QCD vac- 
uum may be studied perturbatively. It may be empha- 
sised that the existence of CEF is itself a consequence of 
the non-perturbative aspects of QCD. Indeed, the chro- 
moelectric field may be considered to be a manifestation 
of the strings in the colour flux tube model [15], which is 
known to provide a natural setting for discussing quark 
confinement [16]. The main difference between this work 
and the conventional approaches is that the latter do not 
take into cognisance the space time dependence of the 
CEF in invoking the production mechanism; they employ 
Schwinger mechanism [17] which is valid only if the field 
is uniform and constant. We make no such assumption 
here (see also [18] in this context). 

Consider the gluons. The lagrangian for pair produc- 
tion may be set up by expanding the gauge potential 
as a sum of the background classical potential and its 
fluctuation (/)^ which is operator valued. Keeping terms 
quadratic in the fluctuations in the Yang -Mills action, 
we get 

+{d^<l>t - 5.<^P(C7''V- + ^''''C-)] + 0{g^). (1) 

Taking into account the Wu Yang ambiguity [19], and 
the studies of Brown and Weissberger [20]. one can ar- 
gue that C° should have an abelian form, at least for 
non vanishing leading order contributions. By a suit- 
able gauge choice, we may write = S/^^ J2i r)Sa,i 
where the summation runs only over the diagonal gener- 
ators of the gauge group. Note that generates only 
an electric component. 

In determining the production rate, as a function of 
time, the crucial step is not to evaluate the S-Matrix. In- 
stead, we study the time dependent evolution of the state 
\^ > {t) in the Fock space using the standard Schwinger- 
Dyson expansion for the Unitary operator U{t,0), with 



the boundary condition \'i> > {t = 0) = \vac >. The in- 
stant t = is singled out as the moment of the creation 
of CEF. The state at any time is then projected on to the 
two gluon sector. The mass shell condition is imposed as 
a constraint on the physical states. 

Denote the two gluon 

state as \gg >= \pi,P2',si,S2',ci,C2 >, labelled by the 
momentum, spin and colour quantum numbers respec- 
tively. In the leading order, the production amplitude 
may be written as 

< 99\Tm >= ^.^-I^^Hn) • 

C'''''{Ei+E2;pi+P2;t) (2) 
where T{t) = U{t,0) - 1. Further, 

(7» = ^*dtie-^(^i+^^)*i J (ffexp{i{p,+p2yf^C°'''ih,r) 

is the incomplete Fourier transform of the gauge fleld 
and Ei are the energies carried by the gluons. The cor- 
responding expression for the qq production is given by 

< qqinm >= J£^-^^CST^\,,^ul{p^)v.,A-P2), (3) 

with \qq >= \pi,p2; si, S2] ci, C2 >. T° are the generators 
of the gauge group in the fundamental representation, 
while u,v are the usual Dirac spinors. 

The probability that a pair is produced any time dur- 
ing the interval (0, i) is given by | < £^\T{t)\vac > p, £, 
standing for either gg or qq. In [1], rates were extracted 
by taking the derivative of the probability with respect 
to time. Here, we study the probabilities directly. 

In order to extract the size parameters, wc deviate from 
the evaluation performed in [1] . We label the two parton 
states by their position coordinates ^1,^2, which entails 
the standard Fourier traansform of Eqns. 2 and 3 in the 
variables pi,p2. The magnitude squared of the ampli- 
tudes now has the significance that it is the probability 
that a pair gets created in the interval {0,t), with the 
particles found at fi,f2 at the instant t. For our pur- 
poses we may sum over the spin, colour and one of the 
position coordinates. The resulting quantity P(r, t) gives 
the (unnormalised) probability density for the parton as 
a function of time. Please note that it is incorrect to 
interpret t as the instant at which the parton is created. 

We employ the same CEF that was introduced in [1] 
to study the prdouction rate in phase space. It is given 
by 

= SMSa,3 + Sa,s)exp{{\z\ - t) /to} 6 {1)6(1'' - z^) 

This choice is close to the boost invariant configuration 
required in the Bjorken scenario [21]. The gauge group 
is SU{3), and the field is dependent only on z,t. The 
probability along the transverse directions is, therefore. 
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uniform, i. e., a gaussian distribution with a width given 
given by the diameter of the nucleus. This dimension 
remains constant in time. However, the longitudinal ex- 
tension has to be determined explicitly. 



III. THE LONGITUDINAL EXTENSION 

We obtained the production probability density for 
gluons and quarks in the one-particle configuration space 
by performing a multi-dimensional fast fourier transform 
on the momentum space amplitude, and integrating over 
the unwanted degrees of freedom. The densities are 
shown in Figs. 1 and 2 as functions of z for different 
time instants namely t = 1, 2, 5, 10, 20 in units of 
tQ. The probability densities arc displayed only on the 
positive z axis. The negative half is symmetric about 
z = 0. The range of the probabilities (unnormalised!) is 
truncated in order to improve the visibility of the impor- 
tant features of the curves. The probability densities at 
z = are large but finite. These curves shed light on an 
interesting aspect of the source. 

Recall that the CEF has a support in the interval 
z G (— i, +t) at any time t, with the two nuclei at \z\ = t 
moving with unit speed in opposite directions. We may 
expect, naively, that the particle production should also 
be restricted within the interval (— The graphs 
clearly bely these classical expectations. The required 
features are most pronounced in Fig. 1, where we show 
the results for the gluons. Strikingly, the probability den- 
sity does not terminate &i z = t. Instead, it extends be- 
yond, as a broad plateau all the way upto z — 2t, where it 
falls abruptly, and almost discontinuously. This feature 
is universal for all the curves, indicating that the point 
z = 2tis naturally chosen as the boundary by the system 
dynamics. In that sense, we do not have to extract any 
length scale by curve fittings. 

The quark results are shown in Fig. 2, again for times 
ranging from i = 1 to < = 20, in units of io- The results 
are not that vivid at earlier times. However, the features 
get more pronounced as the system evolves so that the 
plateau structure is clearly delineated at i = 10/m. 

The remarkable aspect of these conclusions is that the 
size of the system is quintessentially non-classical. In- 
deed, the particles are found in regions where the field 
identically vanishes. Since the boundary of the CEF lo- 
cated at z = ±i is moving with the speed of light, causal- 
ity forbids the particles produced any time in the inter- 
val (0, t) within the range z S (— t) from reaching the 
regions outside the interval. One concludes that the par- 
ticles are indeed produced in the region where CEF van- 
ishes identically! We suspect that this could be a generic 
property of the class of CEF that we have employed. It 
is important to note that the plateau position and its ex- 
tension has a topological nature: the magnitudes change 



smoothly as the sytem evolves with the time, but the lo- 
cation of the shoulder is pegged at the value \z\ = t. A 
straight forward conclusion that one draws is that longi- 
tudinally, the QGP source is twice as large as the space 
between the two nuclei. The fireball extends beyond the 
two nuclei. 
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FIG. 1. Unnormalised longitudinal probability density for 
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FIG. 2. Unnormalized longitudinal probability density for 
quarks as at t = 1, 5, 10 and 20 

A more quantitative understanding can be obtained by 
looking at the relative probability for finding the particles 
in the regions |z| < t and t <\z\ <2t. The ratio for the 
gluons varies from a maximum of about 38% &i t — 1.0, 
to 15% at t = 20. For the quarks, on the other hand, it 
varies from 19% at i = 1.0 to 7% at t = 20. Clearly, the 
non-classical features are the most pronounced when the 
field is varying most rapidly, and gets less important as 
the field evolves in time. Also, should these results seem 
surprising, we may recall that the probabilities in the 
regions where the field is vanishing are, in fact smaller 
than the probability that an oscillator in its ground state 
is found in the classically forbidden region. With this 
perspective, we may conclude that the production in the 
regions where the field vanishes identically is a quantum 



3 



effect; the occurence of the step at \z\ = 2t is perhaps a 
feature of the class of fields which have a boost invariant 
nature, or are close to it. 

A final question remains as to why the spread in the 
quarks is less than that of gluons as indicated by the 
numbers above. The reason for this may be attributed, 
in part, to the fact that the two ghion amplitude is an- 
tisymmetric in the colour and the momentum indices of 
the two gluons (see Eqn. 2. By Bose symmetry, it is 
symmetric in the spins). Consequently, the amplitude 
will be antisymmetric in the position variables ri , r2 as 
well. The quark case on the other hand, has the charge 
conjugation symmetry which, when employed similarly, 
leads to a symmetric behaviour in the spatial part which 
inhibits the spread. In fact, the spread in the momentum 
probabilities is larger for quarks than for the gluons, as 
found in [1]. By the uncertainty principle, one may ex- 
pect the quark size to be smaller than that of the gluons. 

We thank Rajcsh Gopakumar for a discussion and D. 
D. Bhaktavatsala Rao for assistance in preparing the 
manuscript. 
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